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ABSTRACT 

Aims. The aim of this work is to investigate the surface inhomogeneities of a young, late-type star, SAO 51891, at 
different atmospheric levels, from the photosphere to the upper chromosphere, analyzing contemporaneous optical 
high-resolution spectra and broad-band photometry. 

Methods. The full spectral range of FOCESOCAHA (R ~ 40 000) is used to perform the spectral classification and 
to determine the rotational and radial velocities. The lithium abundance is measured to obtain an age estimate. The 
BVRIJHKs photometric bands are used to construct the spectral energy distribution (SED). The variations in the 
observed BV fluxes and effective temperature are used to infer the presence of photospheric spots and observe their 
behavior over time. The chromospheric activity is studied applying the spectral subtraction technique to Ha, Call 
H & K, He, and Ca n IRT lines. 

Results. We find SAO 51891 to be a young K0-1V star with a lithium abundance close to the Pleiades upper envelope, 
confirming its youth (~ 100 Myr), which is also inferred from its kinematical membership to the Local Association. No 
infrared excess is detected from analysis of its SED, limiting the amount of remaining circumstellar dust. We detect 
a rotational modulation of the luminosity, effective temperature, Call H&K, He, and Call IRT total fluxes. A simple 
spot model with two main active regions, about 240 K cooler than the surrounding photosphere, fits very well the 
observed light and temperature curves. The small-amplitude radial velocity variations are also well reproduced by our 
spot model. The anti-correlation of light curves and chromospheric diagnostics indicates chromospheric plages spatially 
associated with the spots. The largest modulation amplitude is observed for the He flux suggesting that this line is very 
sensitive to the presence of chromospheric plages. 

Conclusions. SAO 51891 is a young active star, lacking significant amounts of circumstellar dust or any evidence for 
low mass companions, which displays the typical phenomena produced by magnetic activity. The spots turn out to be 
larger and warmer than those in less active main-sequence stars. If some debris material is still present around the star, 
it will only be detectable by future far-infrared and sub-mm observations (e.g., Herschel or ALMA). The RV variation 
produced by the starspots has an amplitude comparable with those induced by Jupiter-mass planets orbiting close to 
the host star. SAO 51891 is another good example of an active star in which the detection of planets may be hampered 
by the high activity level. 

Key words. Stars: fundamental parameters - Stars: activity - Stars: late-type - Stars: individual: SAO 51891 - Stars: 
planetary systems: protoplanetary disks 



1. Introduction 

In the solar neighborhood there are several stellar kinematic 
groups (SKGs) including stars with a common space mo- 
tion. The origin of these SKGs, like the Local Association 
(LA), or Pleiades moving group, could be the evapora- 
tion of an open cluster, the remnants of a star-forming re- 
gion, or a superposition of several star formation bursts at 
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different epochs in th e adjacent cells of the velocity field 
([Montes et al.l l2001af ). Therefore, the study of late- type 
stars in these young SKGs is important to better under- 
stand the recent star formation history of the solar neigh- 
borhood. 

Moreover, these young solar-mass stars just arrived on 
the Zero- Age Main Sequence (ZAMS), or on their way to 
it, are in a very important phase of their evolution. Indeed, 
as they approach the ZAMS, they spin up while getting 
free from their circumstellar disks, which eventually "con- 
dense" giving rise to planetary systems. At the same time, 
the stars experience angular momentum loss resulting from 
a magnetized stellar wind. In these evolutionary phases, it 
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is therefore very important to define the physical conditions 
that affect the subsequent evolution of the star and any pos- 
sible planets. In addition to basic stellar parameters (T e gp, 
\ogg, [Fe/H]), it is necessary to know the rotation rate and 
the level and behavior of the magnetic activity. 

It is well known that the lower temperature in sunspot 
umbrae and penumbrae as well as in starspots is due to the 
blocking effect of emerging magnetic flux tubes over the 
convective energy transport in the sub-photospheric lay- 
ers. The structure of a magnetic flux tube breaking into 
the stellar atmosphere is traced by the configuration of 
active regions (ARs) at different levels, like photospheric 
spots and chromospheric plages. Therefore, an accurate de- 
termination of spot temperatures and sizes such as that 
derived from the monitoring of line -depth ratios (LDRs) 
along a complete rotation period (jCatalano et alj [2002, 
iFrasca et al.l l2005) and the asso ciated variability of hydro- 
gen, helium and calcium lines (|Frasca et al.l 1200 0. 2008a; 
IBiazzo et al.ll2006l l2007bh can provide crucial information 
on the magnetic field topology in young stars to be com- 
pared with the results for older main sequence (MS) and 
giant stars. 

In the present work we apply this analysis to the 
young star SAO 51891 making use of contemporaneous BV 
photometry and spectroscopic observations. The chromo- 
spheric activity was studied by means of the Ca II H & K 
(A =3968.49 A, 3933.68 A), He (A = 3970.074 A), Ha 
(A =6562.849 A), Can infrared triplet (IRT; A =8498.06 A, 
8542.14 A, 8662.17 A) lines. Our main aim is the determina- 
tion of the spot temperature and filling factor and to infer 
the distribution of chromospheric ARs of SAO 51891. The 
general scope is the comparison of the results obtained for 
very young stars, such as SAO 5 1891, with those found for 
MS and RS CVn st ars (see, e.g.. lFrasca et alj|2005l l2008at 
IBiazzo et al1l2007bh . 

Our work is organized as follows. In Section[2] the main 
stellar parameters of SAO 51891 are briefly discussed. In 
Section[3] we give the details of the observations and data 
reduction. The spectral classification, the rotational and ra- 
dial velocities, the space motion, the spectral energy distri- 
bution, the position on the HR diagram, the lithium abun- 
dance, and the metallicity are reported in Sections 14. 1114. 61 
The diagnostics of photospheric and chromospheric activity 
are analyzed in Sections |4. 71 and l4~51 while the spot model is 
presented in Sect ion [431 Section[5]contains the conclusions. 



2. The case of SAO 51891 (V383 Lac, BD+48 
3686, 1RXS J222007.2+493014) 

SAO 51891 was identified as the optical counterpart of 
the ROSAT extreme ultraviolet source RE J2220+493 
(| Pounds et al.lll99"il iPve et aflll995T). Recen t spectroscopic 
and photometric s tudie s dMulliss fc Boppl 1994 Jeffries 
19951 iHenrv et all Il995l lOsten fc Saarl 119981 : King et al. 
2007) concluded that it is a single active K1V star slightly 
younger than the Pleiades and report ysiai values rang- 
ing from 14 to 20 kms -1 . In particular. IHenrv et alj |l995) 
found the highest peaks in the periodogram of their light- 
curve at P=2.42±0.01 days and 1.70±0.01 days, with the 
former stronger than the latter (both in B and V bands). 
They found an amplitude of 0^08 with significant cycle- 
to-cycle variations of the magnitude at minimum. These 
authors proposed also a minimum radius of 0.77 Rq and an 



inclination i close to 90°. This star does not show any rel- 
evant radial velocity variation as t h e mean value of —20.19 
kms -1 obtained by iMontes et al.l (|2001al ) appears consis- 
tent with the range of values (from —19.4 to —22.1 kms -1 ) 
reported in the literature. This sug gests that SAO 5 1 891 is 
very likely a single star f Table fT|). iMulliss fc Boppl (| 19941) 
found the Ha, Call H and C an IRT lines filled- in with 
emission. Mo ntes et al.l (|2001al) reported notable emission 
in the Ca n H & K and He lines, excess chromospheric emis- 
sion in the hydrogen Balmer lines and emission core rever- 
sal in the Call IRT lines. On one of their observing nights, 
they detected an increase in the excess emission, with broad 
emission wings in the Ha residual profile. They explained 
this variation as probably due to a small-scale flare or to the 
transit of an active region. They also determined a lithium 
equivalent width (EWui) of 257 mA, similar to the values 
of 250 a nd 27 7 mA given by iMulliss k Boppl (|1994| ) and 
Ueffriesl (|1995( ). respectively (Table [1]). This large EW U i 
value, close to the upper envelope of the Pleiades, indicates 
a young a ge. From its kinematic s, SAO 51891 has been con- 
sidered bv lMontes et al.l (|2001al) as a member of the LA. 

Infrared observations place limits on the amount of cir- 
cumstellar dust present (see Section 14.21) . Of course these 
limits only probe dust grain populations with large emit- 
ting areas and do not place constraints on the presence 
of very low mass companions (or planets'). iMetchevI (|2006h 
has conducted an adaptive optics (AOs) survey for faint 
companions in the vicinity of the star. Although several 
candidates were found, no co-moving companions were 
identified, placing limits on brown dwarf companions (> 
50% complete) more massive than approximately the deu- 
terium burning limit at 0.013 Mq beyond approximately 
20 AU. SAO 51891 was also a target of the near-infrared 
AO sear ch for giant planets n amed "Gemini Deep Planet 
Search" (|Lafreniere et a l. 2007). Null results from their sur- 
vey rule out potential planetary mass companions with 
masses greater than 5 Mj beyond 25 AU with 58% proba- 
bility. 

3. Observations and data reduction 

3.1. Spectroscopy 

The spectroscopic observations were conducted in 2006 
at the 2.2-m Cassegrain telescope of the German-Spanish 
Calar Alto Observatory (CAHA, Sierra de Los Filabres, 
Spain) during four nights from 13 to 16 August. The 
Fiber Optics Cass egrain Echelle Spectrograph (FOCES; 
iPfeiffer et ai]|1998D was used with the 2048x2048 CCD de- 
tector Site#ld (pixel size = 24 /im, read-out-noise=5 e~ 
r.m.s., conversion factor = 2.3c~/ADU) allowing us to 
achieve, with an exposure time of 35 minutes, a signal- 
to-noise ratio (S/N) in the range 70-100, depending on 
airmass and sky conditions. The wavelength range covered 
by the 88 echelle orders was 3720-8850 A. We observed 
SAO 51891 at least twice per night obtaining in total 10 
spectra of the object. 

We chose the Site#ld CCD for its higher quantum ef- 
ficiency in the wavelength range of our interest and the 
200 fim fiber-mode in order to avoid substantial light losses 
and maximize the S/N. 

The spectral resolution, determined from the full width 
at half maximum (FWHM) of the emission lines of the Th- 
Ar calibration lamp, was in the range 0.15-0.22 A from the 
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Table 1. 

aturc. 



Physical parameters of SAO 51891 from the liter- 



Table 2. Physical parameters of SAO 51891 derived in this 
work. 



Parameter 



Value 



Reference 



fx (ctss 1 ) 
Spectral Type 

7T (mas) 



H a (masyr x ) 
fis (masyr -1 ) 
vsini (kms -1 ) 



P rot (days) 

R (Re) 
EWui (mA) 



Kad (kms 1 ) 



U Q (kms" 1 ) 
V Q (kms -1 ) 
W Q (kms -1 ) 

T eff (K) 



1.03±0.05 

KOV/IV 

K1V 

42.4±7.2 

36.3±5.0 

20 

93.4±2.5 

5.0±2.5 

15 

16 ± 1 

19.8 

20 

2.42±0.01 

2.410±0.003 

>0.77 

250 

277 

257 

-19.4 

-19.8 ±0.2 

-22.1 ±0.6 

-20.19 ±0.12 

7.06±1.43 

-22.19±0.34 

-3.90±0.86 

5200 



Voges et al. (1999) 
Jeffries (1995) 
Henry et al. (1995) 
Hog et al. (2000) 
Montes et al. (2001b) 
Carpenter et al. (2005 ) 
Hgg et al. (2000) 
Hag et al. (2000) 
Jeffries (1995) 
Henry et al. (1995) 
Fekel (1997) 
Osten fc Saar (1998) 
Henry et al. (1995) 
Xing et al. (2007) 
Henry et al. ( 1995) 
Mulliss fc Bopp (1994) 
Jeffries (1995) 
Montes et al. (2001a) 
Jeffries (1995) 
Henry et al. (1995) 
Osten fc Saar (1998) 
Montes et al. (2001a) 
Montes et al. (2001a) 
Montes et al. (2001a) 
Montes et al. (2001a) 
Osten fc Saar (1998) 



blue to the red, yielding a resolving power R = A/AA ~ 
40 000. 

The data reduction was performed by using the 
echelle task of the IRAlQ package following the stan- 
dard steps of background subtraction, division by a flat- 
field spectrum given by a halogen lamp, wavelength cali- 
bration using the emission lines of the Th-Ar lamp, and 
normalization to the continuum through a polynomial fit. 

We removed the telluric water vapor lines near Ha using 
the spectra of Altair (A7V, vsini = 245 kms -1 ) acquired 
during the observing run, and by applying an interactive 
procedure, which allowed the intensity of the template lines 
to vary (leaving the line ratios unchanged) until a satisfac- 
tory m^t^wjt^ie^h of the observed spectra was achieved 
see iFrasca et al1l2000D . 



3.2. Photometry 

The photometric observations were performed in the B and 
V Johnson filters with the 91-cm Cassegrain telescope at 
the M. G. Fracastoro station (Serra La Nave, Mt. Etna, 
Italy) of the Osservatorio Astrofisico di Catania (OACt). 
The observations were made with a photon-counting re- 
frigerated photometer equipped with an EMI 9893QA/350 
photomultiplier, cooled to — 15°C. The dark current noise 
of the detector, operated at this temperature, is about 1 
count /sec. A 21" diaphragm was used. 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the 
Universities for Research in Astronomy, inc. (AURA) under co- 
operative agreement with the National Science Foundation. 



Parameter 



Value 



Spectral Type 

(B-V) 

V 

Ay 
Prot 

L 

[Fe/H] 

logff SED 

logff SPEC 

rpSEO 

J Gff 

rrnSPEC 

J eff 

T LDR 

J cff 

vsini 
Rsini 

i 

M 

Vrad 

Uq 
Vq 
W @ 
Age 



K0-1±1V 
0.785±0.014 
87M85 ± 0™010 
0.017 

2.62±0.45 days 
0.215 ± 0.075 L e 
0.28±0.05 
4.0±0.1 
4.3±0.1 
5420±150 K 
5350±100 K 
5249±70 K 
19.0±0.5 kms -1 
0.8±0.1P Q 
75°±15° 
O.8±O.2M 
260±10 mA 
-19.67±0.24 kms" 1 
±5.68±1.46 kms" 1 
-21.32±0.46 kms -1 
-3.05±1.01 kms" 1 
< 100 Myr 



SAO 51891 was observed in 2006 from 14 to 21 August 
for a total of 8 nights, using BD±48 3666 (V = 8'. n 308, 
B-V = 0.484) and HD 12028 (V = 8'. n 520, B — V= 0.050) 
as comparison (C) and check (Ck) stars, respectively. 
Additionally, local photometric stand ard stars, as w ell as 
stars of some Landolt's selected areas (lLandolti ri992) were 
also observed during the run in order to determine the zero 
points and the transformation coefficients to the Johnson 
standard system, respectively. The differential magnitudes, 
in the sense of variable minus comparison (V — C), were 
corrected for atmospheric extinction using the seasonal av- 
erage coefficients for the Serra La Nave Observatory. 

A measurement on each object consisted of two inte- 
gration cycles with a 3-sec exposure time in the B and 
V filters, and a typical observing sequence, sky-Ck-C-V- 
C-V, was repeated many times per each night. The data 
were reduced by means of the photometric data reduction 
package PH OT designed for the photo electric photometry of 
the OACt (|Lo Presti fc Marilli|[T993h . The photometric er- 
rors, estimated from measurements of standard stars with 
a brightness comparable to the program stars, are typically 
a v w 0^010 and a B -v « 0^014. 



4. Data analysis and results 

4.1. Spectral classification and rotational velocity 

To derive the spectral type and the vsini of the tar get, we 
used rotfit, a code written by IFrasca et all (|2003l ) in the 
IDL0 environment and successfully applied by us for the 
spectral classification of single-lined ac tive binaries in th e 
RasTyc sample of stellar X-ray sources (jFrasca et al. 2006). 



2 IDL (Interactive Data Language) is a trademark of ITT 
Visual Information Solutions (ITT VIS). 
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The code simultaneously finds the spectral type and the 
v sin i of the target searching for, into a library of standard 
star spectra, the spectrum which gives the best match (min- 
imum of the residuals) with the target one, after the rota- 
tional broadening by convolution with a rotational profile of 
increasing vsm.i at steps of 0.5 kms -1 . During our observ- 
ing run we acquired spectra of only nine standard stars with 
different spectral type. Thus , we decided to use a library of 
ELODIE0 Archive spectra (jPrugniel fc Soubiranl l200lf) of 
87 standard stars well distributed in effective temperature, 
spectral type, and gravity, and in a suitable range of metal- 
licities. For the i>sini determination, we chose as a "non- 
rotating" template the F OCES spect rum of 54Psc (KOV, 
vsini — 2.2 ± 0.8 kms" 1 : [FekelHl997t ) acquired by us dur- 
ing the observing run with the same instrumental set-up. 
This allowed us to avoid small differences between FOCES 
and ELODIE spectra, although the latter have nearly the 
same resolution (R ~ 42 000) as the former. Anyway, the 
average value of vsnii found with ELODIE templates is 
exactly the same as that found with the FOCES spectrum 
of 54Psc (v sini = 19kms _1 ). 

Examples of application of the ROTFIT code to two dif- 
ferent spectral regions are shown in Fig. [TJ where the good 
agreement between observed and standard spectra is evi- 
dent. The spectral type we find is K0-1±1 V and we obtain 
for v smi a best-fit value of 19.0±0.5 kms" 1 (Tabled, near 
to the larger values from the literature. 

4.2. Spectral Energy Distribution 

We defined the spectral energy distribution (SED) in the 
optical/infrared domain with the BVRIJH K S obs erved 
magnitudes after a standard de-reddening |Cox 2000). The 
B and V magnitudes are those obtained by us at OACt 
near the maximum brightness (minimum spot visibility; 
cf. Section [4.7. ip . while the R, I, and JHK S magnitudes 
come from the Na val Observatory Merg ed Astrometric 
Dataset (NOMAD; IZacharias et all 120041) . The Amateur 
Sky Survey fTASS; iDroge et al I l200jl~an d Two Micron 
All Sky Survey (2MASS: ICutri et al.ll2003h catalogues, re- 
spectively. These observed fluxes were compared with those 
obta ined from the low-reso lution synthetic NextGen spec- 
tra (jHauschildt et al]|l999f ) integrated over the pass-bands 
of the filters. The model with T eS = 5400 K, log 5 = 4.0, 
and solar metallicity gives the best fit to the observed SED. 
The radius derived from the fit of the SED and adopt- 
ing the Tycho distance is R = 0.54 ± 0.09 Rq, where the 
error is essentially due to the distance uncertainty. This 
val ue is in very close agr eement with that coming from 
the lBarnes fe Evani (|1976h relation and the Tycho distance 
(R = 0.55 ± 0.19 Rq). However, this value of R is at odds 
with the mass-radius r elation for M S stars. Indeed, accord- 
ing to the tabulation of lCoxl (|2000f) . a K0 V star should have 
a T e g — 5150 K, in agreement with our findings, but a ra- 
dius of 0.85 Rq that is much larger than the value implied 
by the Tycho distance. Moreover, if we combine the avail- 
able determinations of projected rotational velocity, and the 
rotation period of the star, we derive a lower limit for the 
stellar radius, i?sini, of 0.9 Rq (using our determination 
of vsini = 19kms _1 ), or 0.7 Rq (if we adopt the small- 
est value, vsini = 15 kms" 1 from literature). Nevertheless, 



3 This echelle spectrograph of the 1.93-m telescope of the 
Observatoire de Haute- Provence is now decommissioned. 



i f we u se the distance of 27.5 pc reported by iMontes et al.l 
(200l3), the stellar radius turns out to be R = 0.64 i? , 
thus redu cing the discrepancy. T he distance of 50 pc re- 
ported by ICarpenter et alj (|2005l ) would instead increase 
the radius to R rs 1.1 Rq. Hence, we see good reasons to 
suspect that the Tycho distance is underestimated and a 
better distance determination, like that expected from the 
future mission GAIA of ESA, is needed to settle this point. 
In any event, the photospheric temperature is only depen- 
dent on the shape of the SED and the distance has a very 
small effect on it through the interstellar extinction, which 
is, in any case, very low. Interpolating the NextGen fluxes 
at 10 K steps, the best solution in temperature, coming from 
the x 2 minimization, is 5420 K (Tabled]). In Fig. [2J we show 
the observed SED in both linear and log scales, the flux 
computed from the model, and the best-fit synthetic spec- 
trum. For comparison, Carpenter et al. (2008) report an 
effective temperature of 5120 ±130 K with a reddening Ay 
= 0, assuming solar metallicity. 

SAO 51891 was also observed within the Spitzer Legacy 
Progra m "Formation and Evolution of Planetary Systems" 
(FEPS: lMever et al.ll2006h . which is a comprehensive study 
of the evolution of gas and dust surrounding sun-like stars 
from the pre-main sequence (PMS) phase to the age of the 
Sun. Data from 3.6-70 /im plotted in Fig. [2] are from this 
survey. T he ratio of 24 /nm emiss ion to 8 /im emission re- 
ported by ICarpenter et alj (|2008f) is c onsistent with tha t 
expected from the stellar photosphere (jMever et al.| [2008). 
Comparing the upper limit on any pote ntial excess emis- 
sion to the models of iHines et"aLT (|2006l ) suggests a limit 
on the warm dust mass of < 10~ 5 M® in 1-10 /mi grains. 
SAO 51891 is not detected by Spitz er at 70 /im at a sig- 
nificant level (jCarpenter et al.ll2008l ). Comparing the 3(7 
upper limit with the 70 /im flux to the emission expected 
from the star, we can estimate i 71 (70-disk)/F'(70-star) < 16. 
Comparing this limit t o debris disk mode l s of e arly K stars 
(e.g., HD 31392) from iHillenbrand et alj (|2008l ) suggests a 

limit to the cool (50 K) dust mass of < 10 4 Af© in 10 /im 
grains. An upper limit to t he flu x at 850 /im, from a sur- 
vey by iNaiita fe Williams! (|2005l ) using the Submillimiter 
Common-User Bolometric Array (SCUBA) at the James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea (Hawaii), 
is also displayed. More sensitive upper limits on any cold 
dust will require additional far-IR observations with the fu- 
ture mission Herschel of ESA or mm-wave interferometers 
such as ALMA. 



4.3. Position on the HR diagram 

To check the evolutionary status of SAO 51891, we 
placed it onto the Hertzsprung-Russel (HR) diagram 
(Fig. [3]) and compared its p osition with the PMS e volu- 
tionary tracks calculated by iPalla fe Stahlerl |l999) and 
iD'Antona fc Mazzitellil (|1997fh The mean effective temper- 
ature of SAO 51891 was derived considering three determi- 
nations: i) 5420±150K obtained interpolating the NextGen 
flux values at 10 K steps (cf. Section r4~2"|) : ii) the "spectro- 
scopic" value of 5350 ± 100 K obtained imposing the condi- 
tion that the Fel abundance does not depend on the exci- 
tation potentials of the iron lines (cf. Section 14. 6|) ; iii) the 
"LDR-method" value of 5249 ± 70 K corresponding to the 
maximum of the temperature-curve, representative of the 
unspotted (or less spotted) photosphere (the error was eval- 
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Fig. 1. Observed FOCES spectrum of SAO 51891 in the Mgi triplet (upper panel) and 6400 A (lower panel) spectral 
regions together with the ELODIE standard spectrum broadened at 19 kms -1 overplotted with a dashed line. In each 
box the difference (observed — synthetic) is also displayed with a dotted line. 



uated as the quadratic sum of the mean error of « 10 K on 
the individual T c ff values and the rms error of the LDR-T e g 
calibration, cf. Section f4. 7. 2p . At the end, a mean value of 
5350 ± 100 K was adopted. The stellar luminosity was de- 
rived according to the relation L = 4irR 2 aT 2 s , where the 
radius obtained from the SED and assuming the Tycho 
distance is R — 0.54 ± 0.09 i? Q . This yields a value of 
L = 0.215 ± 0.075 L which places SAO 51891 beneath the 
ZAMS (Fig. [3]). Therefore, we calcula ted the lumino s ity also 
al ues rep < 
(|2005T l 

and L = 0.960 L Q , respectively. These values, displayed 
with different symbols in the HR diagram (Fig. [3]), appear 
to be more reliable. We thus conclude that the Tycho dis- 
tance is most probably underestimated. 



usin g the distance va l ues rep orted bv lMontes et alj (|2001bl ) 
and lCarpenter et all (|2005l ) finding L = 0.293 ± 0.084 L Q 



With the latter luminosity values, the position of 
SAO 51891 in the HR diagram is consistent with a ZAMS 
or a post-T Tauri star, in agreement with the lithium con- 
tent and kinematics (cf. Sections 14.41 and 14. 5|) . We derive a 
mass of about 0.9 ± 0.1 M Q . 



4.4. Radial Velocity and space motion 

The heliocentric radial velocity (RV) measurements were 
obtained by means of the cross- correlation technique (e.g., 
lSimkinlll974L iGunn et al.l ll996) using the IRAF task fx- 
COR. A FOCES spectru m of the RV standard star a Ari 
(Whelm = -14.2 kms-^|I^[l967|) was used as a tem- 
plate. In order to take advantage of the wide spectral cov- 
erage offered by FOCES, we cross-correlated each spectral 
order of the SAO 51891 spectra with the template, exclud- 
ing only the orders with low S/N ratio (e.g., the 80th and 
81st, which furthermore include the Can H and K lines) or 
contaminated by broad and/or chromospheric lines (e.g., 
Ha, Nail, and Can IRT) or by prominent telluric features 
(e.g., the O2 series at A ~ 6275 A). Ultimately, about 60 
orders were considered for the cross-correlation function 
(CCF). In order to better evaluate the centroids of the 
CCF peaks, we adopted Gaussian fits. The standard errors 
of the RV values in each order were computed using the 
FXCOR task according to the fitted peak height and the an - 
tisymmetric noise as described by iTonrv fc Davisl (|1979f ). 
For each spectrum, the RV values from individual orders 
were averaged with weights Wi = 1/of (where <Ji is the 
RV error for the z-th order). The resulting RVs are listed 
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Photosphere model: T e {f 



5400 K, logg = 4.00, R = 0.54 R©- 




; A Synthetic fluxes weighted with the BVRIJHK bonds 
O BVRIJHK dereddened fluxes 

5.0X10 3 1.0x10 a 1.5X10'' ; 

A (A) 



□ IRAC/MIPSSSpitzer fluxes 
O SCUBA flux 

— IRS@Spitzer low-resolution spectrun 

- Blockbody flux 




Fig. 2. Top panel: SED of SAO 51891 as deduced from 
BVRIJHKg observed magnitudes (open circles) compared 
to the synthetic fluxes (triangles) by integrating in the 
proper passbands a synthetic NextGen (full line) model 
for T cff = 5400 K and log# SED = 4.0. B and V data are 
referred to our observations, while R, I, and JHK R data 



were taken from N OMAD (IZacharias et alj I2004D TASS 
(jDroge et al.ll2006D . and 2MASS (|Cutri et all 12003 ) cata- 



logues, respectively. Bottom panel: Same figure shown in 
log-log scale, where the IRAC/MIPS @gpiteer (at 3.6, 4.5, 
8, 24, 70 u rn: ICarpenter et al.ll200l and SCUBA/JCMT 
(at 850 /im; iNaiita fc Williama " 20051 ) fluxes are also dis- 
played. The dashed line represents the result obtained con- 
sidering a black body at T B g = 5400 K, the continuous line 
in the 7.6-37 ixm range is instead the lRS@Spitzer spectrum 
(|Carpenter et al.ll2008[) . 



in Table [3] The average RV value over the entire observ- 
ing run is RV = -19.67 ± 0.24 kms" 1 (Tabl e H), which is 
close to the recent determination obtained bv lMontes et al.l 
(|2001af) . As shown in Fig. [5] and in Tabled RV differences 
up to 0.51 kms -1 were measured in our spectra. Since this 
value of RV amplitude is about three times the average er- 
ror on individual measurements, we consider it significant. 
These variations are correlated with the rotational phase 
and are ascribable to line-asymmetry changes caused by 
starspots rather than to an unseen planetary companion 
(see Section g31). 

We used our RV determination combined with the 
parallax and proper motion from the Tycho-2 catalogue 
(Table [I} to derive the Galactic space- velocity components 
(Uq, Vq, Wq) in a left-handed coordinate system (positive 
in the direction of the Galactic anti-center, the Galactic 
rotation, and the North Galactic Pole). We considered 




D'Antona & Mazzitelli (1997) 
Polio & Stohler (1999) 



3.75 
loq Teff 



Fig. 3. HR diagram of SAO 51891. A dot, a diamond, 
and a square are used fo r the values of lum inos- 
ity derived from the Tycho, iMontes et al.l (|2001b[ ). and 
ICarpenter et al.l (120051) distances , respe ctively. The evolu- 
tionary tracks of iPalla &: Stahle"i1 (|l999f ) are shown by con- 
ti nuous lines. Dotted line s are i nstead used for the tracks 
of iD'Antona fc Mazzitellll (|1997fh 



the general outline of iJohnson fc Soderbloml (|1987[ ) with 
the FK4 coordinates at epoch=1950. The uncertainty was 
obtained using the full covariance matrix taking into ac- 
count the error contributions of Vr a d, [is, and tt. The 
values derived are listed in Table [2] and the Boettlinger 
Diagrams in the (U, V) and (W, V) planes are plotted in 
Fig. [4] where the position of some young SKG is also dis- 
pla yed. Our (U,V , W) va lues are very close to those found 
bv iMontes etatl (|2001aD . who derived U = 7.06 ± 1.43 
kms" 1 , V = -22.19 ± 0.34 kms" 1 , and W = -3.90 ± 0.86 
kms -1 as Galactic space- velocity components. We remark 
that the other values of parallax listed in TableQ] do not 
change significantly the values of the space motion compo- 
nents. We als o appl ied a kinematic method developed by 
iKlutsch et alj (|2008l ) to determine the membership proba- 
bility to five young SKGs, namely, the IC2391 supercluster, 
the Pleiades moving group or Local Association, the Castor 
moving group, the Ursa Major group or Sirius superclus- 
ter, and the Hyades supercluster. We find a space velocity 
fully consistent with the young-disk population and a high 
membership probability (64%) to the LA (20-150 Myr). 
Moreover, SAO 51891 is placed inside the subgroup Bl of 
age 20 ± 10 Myr asso c iated with the Pleiades cluster iden- 
tified bv lAsiain etldl (|l999h . 



4.5. Lithium abundance and age 

The lithium abundance was evaluated from the equivalent 
width of the A6707.8 line, EWui- The latter was measured 
using the IRAF task SPLOT, and its error was computed 
by multiplying the integration range and the mean error 
per spectral point evaluated at the continuum on the two 
sides of the line. The lithium abundance was derived us- 
i ng th e curve of growth (COG) method of ISoderblom et ail 
(J 19931 ). The contribution due to the Fei line at A6707.4 
was subtracted using the empirical relationship given by 
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Fig. 4. Space velocities of SAO 51891 in the (U, V) and (W, V) p lanes. The co ntinuous line represents the boundary 
separating young-dis k (YD) and old-disk (OD) stars according to lEggenl (|1996l ). The average velocity components of 
some moving groups (|Eggenlll996l and references therein) are also marked by triangles. 



the same authors: AEW u (rak) = 20(B-V) -3. Lithium 
abundances were t hen corrected for NLT E effects using 
the prescriptions of ICarlsson et all (|l994f ). In the case of 
our target, the LTE lithium abundance turns out to be 
logJV^P «3.18 dex, while the NLTE effects contribute al- 
most halving it (logA^ LTE ss2.95 dex). The lithium abun- 
dance we derive is more than a hundred times the solar 
value (logA£T| = L1 ± .l dex; lLodders et al.ll2009D . 

In the diagram T eS -EWui, SAO 51891 lies above the 
Pleiades upper envelope, indicating an age around one hun- 
dred Myr. This age corresponds to a post-T Tauri (PTT) 
or ZAMS evolutionary stage. Since our target is a sin- 
gle star, its photospheric and chromospheric activity (cf. 
Section | 4.7I) as we ll as its strong coronal X-ray emission 
(jVoges et al.l [1999) should be essentially the effect of its 
young age (see the fx value in Table [T]). 

At this stage, most solar-mass stars are fast or ultra- 
fast rotators, as obser ved in the Pleiades and other yo ung 
open clusters (see, e.g., lMarilli et al.lll997tlBarnesl2003l . and 
references therein), although a fraction of slow rotators is 
also observed. This is attributed to different evolutionary 
processes, such as a larger effectiveness or duration of the 
disk-locking mechanism while the star contracts towards 
the ZAMS. As such, SAO 51891 appears to be somewhat 
similar to the slow rotators in young open clusters. 

4.6. Metallicity 

The iron abundance calculations were performed in LTE 
conditions with an updat ed and impr o ved version of the 
origin al code described in ISpite) (|l967l ). lEdvardsson et al.l 
( 1993?) model atmospheres were used. Equivalent widths of 
spectral lines were measured using the SPLOT task in IRAF. 
The line list and correspon ding atomic data were adopted 
from Pasquini et al. (2004). The effective temperature was 



determined by imposing the condition that the Fe I abun- 
dance does not depend on the excitation potentials of the 
lines. The microturbulence velocity £ was determined by 
imposing that the Fe I abundance is independent of the line 
equivalent widths. The surface gravity log g was determined 
by imposing the Fei/Feii ionization equilibrium. The ini- 
tial value for the effective temperature was the one obtained 
by LDRs (T G L ff DR ; Section 077121) • The initial value of \ogg 
w as obtained taking int o account the evolutionary tracks 
of lPalla fe Stahlerl (| 1999th using the T c L ff DR , solar metallicity 
and M v determined by our photometry and Tycho parallax 
measurements. Due to the good quality of the spectra and 
the many lines present in our wide spectral range, we could 
avoid to use lines in the flat portion of the curve of growth 
(|da Silva et al.|[2~006l ). The initial microturbulence velocity 
was set to 1.0 kms -1 . We derived a metallicity [Fe/H] of 
0.28±0.05, a temperature T^ EC = 5350 ± 100 K, a surface 
gravity log g SPEC = 4.3±0.1 and a microturbulence £ = 1.8 
kms" 1 (Table [J). 

4. 7. Photospheric activity 
4.7.1. Light-curve 

The photometric data acquired contemporaneously to the 
spectroscopic ones allowed us to obtain a light-curve show- 
ing a rotational modulation due to the presence of spots 
(Table IB. 11 Fig. [5]) . The rotational phases were derived ac- 
cording to the following ephemeris: 

HJD^q = 2 453 961.25 + 2 d 42 x E , (1) 

where the initial heliocentric Julian day corresponds to the 
first observing date (namely August 13 st , 2006) at 18:00 UT 
(just before our first observation) and the rota t ional period 
(P rot ) of 2.42 days is taken from lHenrv et all (|1995[ ). This 
period is close to the value of 2.410 ± 0.003 days recently 
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found by IXing et al.l (|2007l ). With a p eriodogram analysis 
and a CLEAN deconvolution algorithm (|Roberts et al.|[T987h 
we find a period of 2.62 ±0.4 5 da ys, which is close t he value 
found bv lHenrv et"aH (|1995l ) and King et"all (|2007ft . taking 
into account the errors. We deci ded to consider the period 
computed bv lHenrv et al.l (|l995l ). because their estimation 
was based on the largest data-set (Table [IJ. 

The V light-curve has a slightly asymmetric shape with 
a peaked top (V max = 8™495) and a variation amplitude 
AF=0™074. The B light-curve has a shape similar to that 
obtained in the V band, with AB = m 079. The B — V 
color shows a scattered modulation with an amplitude of 
only 0'. n 025, which is at a low level of detectability (~ 2a). 
To evaluate the significance of the B — V variation, we per- 
formed a x 2 test by fitting a periodic function ( Fourier poly- 
nomia l) , following the scheme proposed by iBiazzo et alj 
(2006). Notwithstanding the low amplitude, the B—V mod- 
ulation turns out to be significant. However, it can hardly 
be used to constrain the spot parameters and an additional 
information, such as that provided by the line-depth ratios, 
must be used for this purpose. 



It 




4.7.2. Effective temperature from line-depth ratios 

It has been demonstrated that line-depth ratios (LDRs) are 
powerful tools for detecting temperature rotational mod- 
ulations in stars with moderate magnetic activity (e.g., 
Biaz zo et al.l l2007bT) and with a high level of a ctivity (e.g. 



Catalano et al.ll2002t IFrasca et al.1 12005. 20083). Such diag- 



nostics allow us to detect temperature variations as small 
as 10 K at the resolution of o ur spectra and with a sign al-to- 
noise ratio higher than 100 (jGrav fc Johansonlll99lh . The 
precision of this method is improved by averaging the re- 
sults from several line pairs. In particular, for SAO 51891 
we used fifteen line pairs in the optical spectral range 
6200-6300 A and th e LDR-T e g- calibrations developed by 
IBiazzo et all (|2007al) at vsmi = 20 kms . Thus, we ob- 
tained the temperature-curve as the weighted average of 
all the temperature-curves from each LDR (Fig. [5]). The 
standard error of the weighted average was computed on 
the basis of the errors in each LDR-derived temperature. 
As shown in Fig. [5l the average effective temperature is 
in phase with the light-curve, with only a possible shift 
of the maximum by less than 0?1. However, we have a 
gap in the decreasing part of the temperature-curve which 
casts doubts on the reality of this shift. The flat tempera- 
ture minimum, with a dip around phase 0?0, coincides in 
phase with the minimum in the light-curve. This fair coher- 
ence confirms the hypothesis of cool spots as the primary 
cause of the observed variations. The temperature varies 
between 5162 K and 5249 K, with a full amplitude of about 
90 K, which is intermediate betw een the values of ps40K 



found in moderately active stars (IBiazzo et al. 2007rj) and 
ps 130 K found in highly a ctive stars ( Catalano et al.ll2002l : 
iFrasca et al1l2005l l2008ah . 



4.8. Chromospheric activity 

The wide wavelength range of FOCES spectra allowed us to 
study the chromosphere of SAO 51891 by using several lines 
from the near UV to the NIR wavelengths (namely, Call 
H & K, Ha, Can IRT), which carry information on different 
atmospheric levels, from the region of temperature mini- 



Fig. 5. Observed (dots) and synthetic (dashed lines) V, B, 
T e g, Vrad curves of SAO 51891 versus the rotational phase. 



mum to the upper chromosphere. To derive the chromo- 
spheric losses, we used the "spectral synthesis" technique, 
based on the comparison between the target spectrum and 
an observed spectrum of a non-active standard star (called 
"reference spectrum" ) . The difference between the observed 
and the reference spectrum provides, as residual, the net 
chromospheric line emission, which can be integr ated to 
find t he total radiative losses in the line (see. e.g..lHerbig 
1 19851; iBardenl [19851 IFrasca fc Cat aland [l99l iMontes et at 
11995ft . 

The non-active star used as a reference for the spec- 
tral subtraction is 54 Psc (=HD 3651), a K0V star [B - 
V = 0.85) with a very low activity level, a s indicated by 
the l ow value of the Can S index (0.195; iDuncan et alj 
119911 ) and by the low flux in the CaiiK line me asured on 
high- resolution spectra (namely logTrFx = 4.98: ICataland 
1979). This star was also observed during the same run as 
SAO 51891. In Fig. [BJ we show an example of a spectrum 
of SAO 51891 in the Ha, and Can IRT regions, together 
with the standard-star spectrum rotationally broadened to 
vsini = 19 kms -1 which mimics the active star in absence 
of chromospheric activity. The Can H&K region is dis- 
played in Fig. [71 The Ha and Can IRT profiles are clearly 
fillcd-in by emission, with the calcium lines displaying an 
emission core. The Ca n H & K cores exhibit strong emission 
features and He emission is also clearly visible (Fig. [7J. 

The residual equivalent width (EW) of the lines has 
been measured by integrating all the emission profile in 
the difference spectrum (see bottom of each panel of Fig. [Fj] 
and bottom panels of Fig. [7]). The error on EW {(Jew) 
was evaluated by multiplying the integration range by the 
photometric error on each point. The latter was estimated 
by the standard deviation of the observed fluxes on the 
difference spectra in two spectral regions near the line. 
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Fig. 7. Top panel: Example of an observed, continuum-normalized spectrum of SAO 51891 (solid line) in the Can H & K 
region together with the non-active stellar template (dotted line). Bottom panels: The difference spectrum in the Can K 
and H regions. He emission is clearly visible both in the raw and in the difference spectrum. 



4.8.1. Ha line 

It has been widely shown, both from theoretical and ob- 
servational points of view, that the Ha line is one of the 
most useful and easily accessible indicators of chromo- 
spheric emission related to solar and stellar activity. Being 
its source function photoionization- dominated in the quiet 
chromospheres of the Sun and solar-like stars, this line con- 
tains valuable informations on the chromospheric structure 
and it is very sensitive to the strong non-thermal velocity. 
The Ha absorption in the spectra of cool stars provides 
evidence for the existence in these stars of chromospheres 
with significant optical depth in that line. As a consequence, 
Ha is useful for detecting chromospheric solar and stellar 
plages, thanks to their high contrast against the surround- 
ing chromosphere. 

The Ha line of SAO 51891 is always in absorption with 
a strong filling-in of the core that appears broader than 
that observed in MS stars with mild activity like K\ Ce t 
(|Biazzo et all l2007bl) or HD 206860 (jFrasca et all 120001) . 
Moreover, the residual Ha profile of SAO 51891 is almost 
symmetric and does not show a significant variation (see 
Fig. [5] and Appendix [5} ■ 



The values of the residual emission EWn a with their 
errors are listed in Table [3] and plotted as a function 
of the rotational phase in Fig. [8] No rotational modula- 
tion of EWhq, emerges over the data scatter. This result 
maybe is surprising, because rotational modulation of the 
Ha emission has been frequently detected in several highly- 
and mildly-active stars (see, e.g . , Frasca et al.ll2000L |2008bt 
lAlekseev k Kozlova|[2002l [2003; Bi azzo et alj 2006. 2007bh. 



However, in some cases, no clear rotational modulation of 
chromospheric diagnostics has been detected despite the 
contemp oraneous wave-like be haviour in the photometry 
(see, e.g.. lCatalano et al"1l2000D . In such cases, the distribu- 
tion of the active regions, as revealed by the Ha emission, 
may be more homogeneously distributed in the base of the 
chromosphere than in the photosp here or in other atmo- 
spheric layers, as proposed, e.g., by iGarcia- Alvarez! {2003). 
In addition, the Ha line may originate in other physical 
processes, like flares and prominences. 

Since we observe a scatter in the Ha EW of SAO 51891 
larger than the typical data errors, we suppose that addi- 
tional phenomena, such as microflaring, producing intrinsic 
variations of the emission can affect the line profile, hiding 



10 



K. Biazzo et al.: The young active star SAO 51891 




Table 4. Radiative chrome-spheric losses. 



Line 


Flux 




(erg cm -2 s _1 ) 


Ha 


3.0xl0 6 


Can H+K 


4.7xl0 6 


Can IRT 


4.9xl0 6 


Mg II h+k 


4.9xl0 6 



X (A) 

Call AA8498.8542 




8520 8530 
X (A) 

Call A8662 




8655 8660 8665 8670 

X (A) 

Fig. 6. Top of each panel: Examples of observed, 
continuum- normalized spectra of SAO 51891 (solid line) in 
the Ha and Ca n IRT regions together with the non-active 
stellar template (dotted line). Bottom of each panel: The 
difference spectra of the two upper spectra. 



any low-amplitude rotational modulation. The rotational 
modulation of Ca n H&K and He line fluxes (cf . 14.8.21 and 
I4.8.3P supports this idea. 

We also calculated the chromospheric radiative losses 
in the Ha line, fket < following the guidelines by 
iFrasca fe Catalanq (|1994h . i.e. by multiplying the net 
equivalent width, EWn a , by the continuum surface flux at 
A = 6563 A. The latter was evaluated for all the stars of the 
spectrophotometric atlas of iGunn fe Strvkerl (fl98l with 
the an gular diameters calculated from the iBarnes fe Evans! 
(|1976[ ) relation. The continuum flux of SAO 51891 was 
found by interpolating the values for the IGunn fe Strvkerl 
(|1983f) stars at B-V = 0.785. The line flux F Ha is reported 
in Table [4] together with those of other chromospheric di- 
agnostics (see Sections 14.8.21 and 14.8.31) . 



4.8.2. Can H&K and He lines 

The Ca II H & K lines show very strong emission cores with- 
out any detectable self-absorption. They also appear very 
symmetric in all spectra. The He emission is clearly visible 
at the right side of the Can H line, both in the observed 
and in the residual spectra (Fig. [7]) . 

Since the absorption wings of each of the two lines span 
two echelle orders, it was necessary to merge them before to 
proceed with the spectral subtraction analysis. For the or- 
der merging and normali zation we used the same procedure 
as in lFrasca et al.l ()2000f ) . The spectrum of the low-activity 
star 54 Psc, broadened at vsin«=19kms~ 1 , was used as a 
non-active template, as explained in Section [4.81 Although 
a tiny emission/filling- in of the core of the K line is barely 
visible, it is negligible compared with the huge emission of 
SAO 51891. Moreover, we preferred to use the same tem- 
plate for all the activity diagnostics. 

We measured the equivalent widths by integrating the 
emission profiles in the subtracted spectra, as for Ha. For 
deblending the Can H and He lines, we performed fits of 
two Voigt functions by means of the SPLOT task of IRAF. 
Both the H and K EWs display a fair rotational modula- 
tion, nearly anti-correlated with the light-curve. In Fig. [5J 
we show the sum of the net equivalent widths, EWq^. A 
more pronounced modulation is observed in the net equiv- 
alent width of the He line. 

iTuroval (|1994l ) observed a long-lasting He emission (ss 6 
days) above umbral regions for a sunspot group. This is 
a very rare phenomenon in the Sun, observed sometimes 
during flares. Based on the intensity ratio He/Can H, she 
suggested that the sunspot group behaved like a dKe or 
dMe star, and concluded that the emission in the He 
line is related to an enhancement of collisional transitions 
owed to a very large t e mperature gradi ent in the chromo- 
sphere (lFosburv|[l974D . IFosburvl (|1974f) stressed the fact 
that the emission in He could be due to a photoionization- 
recombination driven by the Ca n H line and background 
continuum radiation for stars with low chromospheric den- 
sity, like in Arcturus, but the line is collision-dominated in 
dKe and dMe where the electron density is much higher. 
Moreover, when the electron density in the chromosphere is 
high, the Ca n emission cores are saturated and the relative 
intensity In/ Ik approaches unity. In these conditions, the 
He emission beco mes more prom inent compared to that of 
Can H. Fig. 2 of IFosburvl (fl97l displays the Can H&K 
region of AD Leo, in which the He emission peak is more 
than one third of the Can H intensity (« 0.35 - 0.40). For 
SAO 51891 we find values in the range 0.16 — 0.19 with a 
possible (very scattered) modulation. Hence, if the chro- 
mosphere above an active region is strongly non-thcrmally 
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HJD 


Phase 


v rad 








(A) 


EW He 


(+2 400 000) 




(kms" 1 ) 


(K) 


(A) 


(A) 


(A) 


53 961.379 


0.053 


-19.60±0.19 


5174±12 


0.579±0.036 


1.424±0.055 






53 961.473 


0.092 


-19.63±0.17 


5162±7 


0.633±0.028 


1.394±0.054 


1.168±0.059 


0.175±0.023 


53 961.566 


0.131 


-19.78±0.18 


5188±14 


0.615±0.038 


1.406±0.051 


1.195±0.055 


0.226±0.022 


53 962.453 


0.497 


-19.80±0.17 


5249±10 


0.625±0.030 


1.404±0.054 


1.116±0.046 


0.141±0.015 


53 962.609 


0.562 


-19.57±0.18 


5242±8 


0.572±0.036 


1.409±0.059 


1.036±0.068 


0.116±0.019 


53 963.402 


0.889 


-19.33±0.17 


5204±9 


0.603±0.023 


1.361±0.053 


1.107±0.052 


0.134±0.017 


53 963.555 


0.952 


-19.65±0.19 


5167±13 


0.656±0.036 


1.337±0.067 


1.228±0.069 


0.143±0.019 


53 964.371 


0.290 


-19.72±0.21 


5177±2 


0.560±0.035 


1.431±0.071 


1.289±0.119 


0.205±0.042 


53 964.477 


0.333 


-19.70±0.17 


5201±13 


0.594±0.021 


1.419±0.060 


1.209±0.060 


0.144±0.015 


53 964.574 


0.374 


-19.84±0.18 


5220±9 


0.673±0.028 


1.464±0.060 


1.230±0.053 


0.169±0.014 



heated and the density is high enough, He goes into emis- 
sion and enhances much more than the Can lines do. This 
could explain the high amplitude of the He modulation and 
suggests this line as a sensitive marker of chromospheric 
surface features. 

We also calculated the radiative chromospheric losses 
in the H&K lines analogously as done for Ha. In partic- 
ular, we used two 10 A- wide bands centered at 3910 and 
4010 A, i.e. at the two sides of the H&K lines to per- 
form the flux cali bration, following the prescriptions by 
iFrasca et all (|2000[ ). 

Finally, we evaluated the radiative losses in the h & k 
lines of Mgn on a low-resolution UV spectrum, namely 
LWP26439LL.FITS, the only one available in the IUE0 
Final Archive, and using as non-active template a resam- 
pled IUE spectrum of 54Psc. The radiative losses in the 
Ca II H & K and Mg II h & k, which are also listed in Tablegl 
turn out to be nearly the same. 

4.8.3. Can IRT lines 

These lines share the upper level of the H and K transitions 
and are useful for stellar chromospheric activity studies. 
Their extended wings probe a wide range of photospheric 
layers and are sensitive to the temperature distribution in 
the atmosphere of the star, while their cores are formed 
in the uppermost atmospheric layers and are sensitive to 
the physical conditions in the chromosphere. Some investi- 
gations based on empirical chromospheric models applied 
to active stars on a wide range of activity levels suggest 
that the contribution of the Call IRT lines can be up to 
twice larger than the contribution of the Ca n H & K lines 
(|Dempsev et al.llT993D . As in the Sun, the calcium lines sig- 
nificantly contribute to the total chromospheric losses of 
the active stars, providing useful in formation on th e energy 
balance of stellar chromospheres (|Busa et al.ll2007l . and ref- 
erences therein). 

The lines of the Can infrared triplet present some ad- 
vantages compared to the Ca n H & K lines. First of all, the 
triplet lines lie in a spectral region with a well-defined con- 
tinuum, making the normalization easier during the data 
reduction. Moreover, they are not significantly affected by 
telluric lines and are less affected by atmospheric extinc- 
tion than the visible and ultraviolet lines. Finally, the high 
sensitivity of the new detectors to the near-IR makes these 
lines more easily observable than in the past. For this rea- 
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son the calcium triplet has become a very versatile tool and 
the spectral region 8480-8740 A has been selected for the 
medium resolu tion spectrograph of the GAIA mission, as 
pointed out bv lBusa et"aT1 (|2007l ). 

In SAO 51891 the Call IRT lines are always filled-in 
with a central emission peak never reaching the local con- 
tinuum. The profiles of the two lines AA 8542.14,8662.17 
are almost symmetric, while the A8498.06 line displays an 
asymmetric profile fFigs. lA~2| IA.3jl . 

The net equivalent width of these three lines (in partic- 
ular the A8498.06 line) shows a detectable rotational modu- 
lation, which becomes more evident if we consider the total 
emission of the triplet (see in the bottom panel in 

Fig. [H]). To evaluate the significance of the BWj^ varia- 
tion, we performed the x 2 test, as done for the B — V index. 
The reduced x 2 of the fit is 0.19 for the Fourier polynomial 
fit and 0.34 for the constant function. The probability that 
a fit of 5 free parameters gives x 2 < 0.19 is 0.08%, while 
with one degree of freedom we obtain 44%. This makes the 
rotational modulation of EWq^ significant. The chromo- 
spheric radiative losses in the three lines of the Ca n IRT 
(Table 2]) have been calculated as for the Ha and Call 
H feK lines. They are n early the same as those reported 
by iMontes et alj (|2001af l and turn out to be only slightly 
larger than those we found for the H&K lines. 

4.9. Modeling the light and temperature curves 

In IFrasca et all ([2005) we showed that, with an IDL spot 
model applied to contemporaneous light and temperature 
curves, it is possible to reconstruct the distribution of 
starspots and to remove the degeneracy of the spot param- 
eters temperature and area. We modeled the observed light 
and temperature curves by assuming two active longitudes 
on the stellar photosphere sketched by two circular regions 
whose flux contrast (F sp /Fph) can be evaluated thr ough the 
Planc k spectral energy distribu tion, the ATLAS9 dKuruczl 
H99l and PHOENIX NextGen (lHauschildt et al.lll999D at- 
mosphere models. Even if the real shape of the active re- 
gions can be very different, this approximation is useful to 
define their main average p arameters, like posi tion, relative 
area, and temperature. In IFrasca et alj (|2005l ) we demon- 
strated that both the atmospheric models (ATLAS9 and 
NextGen) provide values of the spot temperature (T sp ) and 
area coverage (A IC \) in close agreement, while the black- 
body assumption for the SED leads to underestimate the 
spot temperature. Since we have no long-term record of the 
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Table 5. Parameters of the two photospheric active re- 

, , , , , , , , , , , , , , , gions. 




8.60 - 



0.0 0.5 1.0 

Phase 

Fig. 8. From top to bottom. EW Ha , EW^ m EW^+f, 
EWii e , and V magnitude as a function of the rotational 
phase. 



photospheric temperature, we assumed the maximum value 
obtained during our observations as the temperature of the 
unspotted photosphere. 

The application of our spot model requires the knowl- 
edge of the inclination of the rotation axis with respect to 
the line of sight. From our SED and the Tycho distance, 
we estimated a radius R = 0.54 ± 0.09 R Q (Table©. By 
combining it with the rotation period of 2.42 days and the 
vsini of 19kms determined by us, we obtain a sini > 1, 
testifying an inconsistency between these parameters. We 
think that the small radius value, coming from the pos- 
sibly underestimated Tycho distance, is the source of the 
inconsistency, because the setback remains also with the 
lowest value of v sinz ~ 15 kms -1 from the literature. If we 
use a typi cal radius for a K0V ZAMS star (R = 0.85 R Q ; 
ICoxl [20001 we obtain an inclination i of nearly 90°, while 
adopting our highest radius esti mate, R = 1.1 Rp>. based on 
the distance of 50 pc reported bv lCarpenter et all ({20081 ) . a 
value of i « 60° is found (Table [J). Thus, we adopted three 
values of i (namely 60°, 75°, 90°) for our spot model. Then, 
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from our unspotted temperature and magnitude, namely 
T ph = 5249 K and V = 8™502 we obtained, for each i, a 
grid of solutions for the temperature curve and another 
one for the light curve. The limb-darkening coefficients 
M6200 = 0-55 and /zy = 0.74 were u sed for the temperature 
and t he V light curve, respectively (| Al-Naimivlll978l : [Claretl 
2000). The intersection of the T e s and V grids provides 
us with the values of spot temperature and area. Figure [9] 
shows that, for any value of spot temperature (T sp /T p h), the 
size of the active regions must increase with the inclination 
in order to reproduce the observed variations. The effect is 
more pronounced for the temperature grid, so that we can 
have an intersection (a solution) only for i = 60°, which 
occurs for T sp /T p h = 0.955 and A IC \ = 0.149. This seems to 
support a low inclination for the rotation axis, i.e. a large 
stellar radius. However, due to the strong uncertainty in the 
inclination, we are only able to give a very rough estimate 
of the spot temperature and area, while the longitudes of 
the two active regions (Table©, corresponding to 0.97 and 
0.29 phases, are reliable. 

We have also applie d a plage model written in IDL (see, 
e.g.. lFrasca et al.l l2000) to the modulation of the total Call 
IRT and HK emission with the aim of gaining some infor- 
mation about the surface inhomogeneities at chromosphcric 
level. Given the scatter in the data, our code allowed us 
only to assess that the BW^f and EW^ 1 curves are 
fairly reproducible with a single plage (A IC \ ~ 0.055) placed 
between the two spots and closer to the smaller one. The 
only meaningful parameter that can be deduced from these 
data is the longitude of the chromospheric AR (w 110°). 
However, this does not exclude two plages around the pho- 
tospheric spots. 

Our spot model allows us also to calculate the radial 
velocity variations induced by the line (and CCF) distor- 
tions owing to the passage of spots (or pla ges) over the 
star di sk if the v sin i of the star is known (see lFrasca et al.l 
2008b]). With the spot parameters listed in Table we cal- 
culated the expected RV curve and we superimposed it to 
the observed RV variations (with a full amplitude of about 
500 ms _1 ) in the bottom panel of Fig. [51 The agreement 
with the data is apparent. This means that it is not neces- 
sary to invoke a low-mass companion (brown dwarf or giant 
planet in a close orbit) to explain the observed variations. 
Moreover, such variations prevent us to detect a Jovian 
planet around SAO 51891, because the RV curve due to 
the planet with an amplitude of less than 100 mr 1 would 
be lost inside the RV variations produced by the starspot. 
Thus, very active stars, such as SAO 51891, are not suitable 
targets for exo-planet search with the RV technique, unless 
a simultaneous detailed study of the magnetic activity phe- 
nomena is carried out. 
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Fig. 9. Grids of solutions for SAO 51891 obtained for three different values of i. The filled circles represent the values of 
spot temperature and area for the light-curve solutions, while the diamonds represent the solutions for the temperature- 
curve. The green lines in the left panel delimit the locus of the allowed solutions accounting for the data uncertainty. 



5. Conclusion 

In this paper we analyzed high-resolution echelle spectro- 
scopic observations and contemporaneous BV photometry 
of the young star SAO 51891. We report an updated spec- 
tral type, revised astrophysical parameters such as T e ff, 
log g, [Fe/H], rotational and heliocentric radial velocity, 
space motion, spectral energy distribution, and lithium 
abundance. Our main goal was to investigate short-term 
variability due to magnetic activity at both photospheric 
and chromospheric levels. 

The kinematics, the lithium abundance, and the level 
of photospheric and chromospheric activity indicate an age 
similar to that of the Pleiades cluster. We also confirm the 
membership in the Local Association, with an age of the or- 
der of 100 Myr. The spectral energy distribution provides 
evidence for lacking of significant amounts of circumstel- 
lar dust and upper limits are derived. The small-amplitude 
(« 500 ms _1 ) radial velocity variations measured from our 
spectra are correlated with the stellar rotation and fully ex- 
plained with the line distortions produced by the starspots. 
Thus there is no evidence for a stellar or sub-stellar com- 
panion, based on radial velocity measurements or direct 
imaging. Our analysis demonstrates the difficulty of detect- 
ing sub-stellar companions (brown dwarfs or giant planets) 
around very active stars from radial velocity measurements. 

From our spectra, covering one and a half stellar ro- 
tation, we find conspicuous chromospheric emission in the 
Can H&K, Can IRT and He lines. The cores of the Ha 
line are also clearly filled in by chromospheric emission. 
All these features confirm the strong magnetic activity of 
SAO 51891. 

We detected a clear modulation of the V and T g curves 
due to spots, and a low-amplitude modulation of the Call 
IRT, H&K, and He emissions ascribed to chromospheric 



inhomogencities. Astonishingly, we did not find any clear 
modulation in Ha, in contrast to what we observed for sev- 
eral other active stars. We speculate that the modulation 
produced by plages is possibly hidden by other activity sig- 
natures, like microflares, which can significantly affect the 
Ha emission. 

The simultaneous solution of the V light-curve and the 
temperature modulation with our spot model allowed us to 
reconstruct the approximate spot distribution on the stellar 
photosphere by adopting an inclination i = 60°. We found 
two fairly large active longitudes with a temperature close 
to the "unspotted" effective temperature (AT ~240K), 
different to what we have found for other mildly-active main 
sequence stars. We want to stress that our model gives us 
only a "rough" estimation of the inhomogeneities parame- 
ters mainly for the following reasons: i) there is a phase gap 
in our data (in particular the spectroscopic observations) 
just after the maximum of the curves which can introduce 
some uncertainty on the "unspotted" level of temperature; 
ii) uncertainties on the distance translate into uncertainties 
in the stellar radius and, consequently, the inclination angle 
cannot be constrained with sufficient accuracy. 

To date, photometric and spectroscopic analysis com- 
parable to that reported here has been conducted on only 
a handful of young stars. In order to place our own so- 
lar system in context, understand the range of planet for- 
mation outcomes as a function of stellar parameters, as 
well as investigate the early angular momentum evolution 
of sun- like stars, more observations of this kind are re- 
quired. For this reason, we already started a program of 
high-resolution FOCES@CAHA and SARG@TNG echelle 
spectroscopic observations of young late-type stars to im- 
prove our knowledge of photospheric/chromospheric inho- 
mogeneities of this stellar population. 
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Appendix A: Spectra of SAO 51891 in the Ha and 
Can IRT regions, and calcium EWs values. 

Ha 



Observed ( ) and synthetic (. . ) spectra Difference spectra 




6549 6558 6567 6576 6549 6558 6567 6576 

A (A) A (A) 



Fig. A.l. Left panels: observed and continuum- normalized 
spectra of SAO 51891 in the Ha region together with the 
non-active template. Right panels: difference between ob- 
served and template spectra. 
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HJD Phase 






EW^f 2 






(+2 400 000) 


(A) 


(A) 


(A) 


(A) 


(A) 



53 961.379 
53 961.473 
53 961.566 
53 962.453 
53 962.609 
53 963.402 
53 963.555 
53 964.371 
53 964.477 
53 964.574 



0.053 
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0.131 
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0.562 
0.889 
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0.290 
0.333 
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Table B.l. V and B — V values. 
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B - V 


(+2 400 000) 


(mag) 
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Fig. A. 3. Left panels: observed and continuum-normalized 
spectra of SAO 51891 in the Ca II IRT A8662 region together 
with the non-active template. Right panels: difference be- 
tween observed and template spectra. 
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